Building pyramids against the evolutionary emergence of pathogens.

Mutations allowing pathogens to escape host immunity promote the spread of infectious diseases in heterogeneous host populations and can lead to major epidemics. Understanding the conditions that slow down this evolution is key for the development of durable control strategies against pathogens. Here, we use theory and experiments to compare the efficacy of three strategies for the deployment of resistance: (i) a mixing strategy where the host population contains two single-resistant genotypes, (ii) a pyramiding strategy where the host carries a double-resistant genotype, (iii) a combining strategy where the host population is a mix of a single-resistant genotype and a double-resistant genotype. First, we use evolutionary epidemiology theory to clarify the interplay between demographic stochasticity and evolutionary dynamics to show that the pyramiding strategy always yields lower probability of evolutionary emergence. Second, we test experimentally these predictions with the introduction of bacteriophages into bacterial populations where we manipulated the diversity and the depth of immunity using a Clustered Regularly Interspaced Short Palindromic Repeats-CRISPR associated (CRISPR-Cas) system. These biological assays confirm that pyramiding multiple defences into the same host genotype and avoiding combination with single-defence genotypes is a robust way to reduce pathogen evolutionary emergence. The experimental validation of these theoretical recommendations has practical implications in various areas, including for the optimal deployment of resistance varieties in agriculture and for the design of durable vaccination strategies.

Complete genome of Escherichia phages Carena and JoYop.

Escherichia phages Carena and JoYop were isolated from water samples in Abidjan (Cote d'Ivoire). Their genomes comprise 39,283 and 169,193 bp, encoding 44 and 246 predicted genes, respectively. Carena shares 93.4% nucleotide identity with Escherichia podophage CarlSpitteler (Berlinvirus), and JoYop shows 95.6% identity with Escherichia myophage ADUt (Tequatrovirus).

Longitudinal Study of Lactococcus Phages in a Canadian Cheese Factory.

The presence of virulent phages is closely monitored during cheese manufacturing, as these bacterial viruses can significantly slow down the milk fermentation process and lead to low-quality cheeses. From 2001 to 2020, whey samples from cheddar cheese production in a Canadian factory were monitored for the presence of virulent phages capable of infecting proprietary strains of Lactococcus cremoris and Lactococcus lactis used in starter cultures. Phages were successfully isolated from 932 whey samples using standard plaque assays and several industrial Lactococcus strains as hosts. A multiplex PCR assay assigned 97% of these phage isolates to the Skunavirus genus, 2% to the P335 group, and 1% to the Ceduovirus genus. DNA restriction profiles and a multilocus sequence typing (MLST) scheme distinguished at least 241 unique lactococcal phages from these isolates. While most phages were isolated only once, 93 of them (out of 241, 39%) were isolated multiple times. Phage GL7 was isolated 132 times from 2006 to 2020, demonstrating that phages can persist in a cheese factory for long periods of time. Phylogenetic analysis of MLST sequences showed that phages could be clustered based on their bacterial hosts rather than their year of isolation. Host range analysis showed that Skunavirus phages exhibited a very narrow host range, whereas some Ceduovirus and P335 phages had a broader host range. Overall, the host range information was useful in improving the starter culture rotation by identifying phage-unrelated strains and helped mitigating the risk of fermentation failure due to virulent phages. IMPORTANCE Although lactococcal phages have been observed in cheese production settings for almost a century, few longitudinal studies have been performed. This 20-year study describes the close monitoring of dairy lactococcal phages in a cheddar cheese factory. Routine monitoring was conducted by factory staff, and when whey samples were found to inhibit industrial starter cultures under laboratory conditions, they were sent to an academic research laboratory for phage isolation and characterization. This led to a collection of at least 241 unique lactococcal phages, which were characterized through PCR typing and MLST profiling. Phages of the Skunavirus genus were by far the most dominant. Most phages lysed a small subset of the Lactococcus strains. These findings guided the industrial partner in adapting the starter culture schedule by using phage-unrelated strains in starter cultures and removing some strains from the starter rotation. This phage control strategy could be adapted for other large-scale bacterial fermentation processes.

In through the Out Door: A Functional Virulence Factor Secretion System Is Necessary for Phage Infection in Ralstonia solanacearum.

Bacteriophages put intense selective pressure on microbes, which must evolve diverse resistance mechanisms to survive continuous phage attacks. We used a library of spontaneous Bacteriophage Insensitive Mutants (BIMs) to learn how the plant pathogen Ralstonia solanacearum resists the virulent lytic podophage phiAP1. Phenotypic and genetic characterization of many BIMs suggested that the R. solanacearum Type II Secretion System (T2SS) plays a key role in phiAP1 infection. Using precision engineered mutations that permit T2SS assembly but either inactivate the T2SS GspE ATPase or sterically block the secretion portal, we demonstrated that phiAP1 needs a functional T2SS to infect R. solanacearum. This distinction between the static presence of T2SS components, which is necessary but not sufficient for phage sensitivity, and the energized and functional T2SS, which is sufficient, implies that binding interactions alone cannot explain the role of the T2SS in phiAP1 infection. Rather, our results imply that some aspect of the resetting of the T2SS, such as disassembly of the pseudopilus, is required. Because R. solanacearum secretes multiple virulence factors via the T2SS, acquiring resistance to phiAP1 also dramatically reduced R. solanacearum virulence on tomato plants. This acute fitness trade-off suggests this group of phages may be a sustainable control strategy for an important crop disease. IMPORTANCE Ralstonia solanacearum is a destructive plant pathogen that causes lethal bacterial wilt disease in hundreds of diverse plant hosts, including many economically important crops. Phages that kill R. solanacearum could offer effective and environmentally friendly wilt disease control, but only if the bacterium cannot easily evolve resistance. Encouragingly, most R. solanacearum mutants resistant to the virulent lytic phage phiAP1 no longer secreted multiple virulence factors and had much reduced fitness and virulence on tomato plants. Further analysis revealed that phage phiAP1 needs a functional type II secretion system to infect R. solanacearum, suggesting this podophage uses a novel infection mechanism.

Genome Sequence of SN1, a Bacteriophage That Infects Sphaerotilus natans and Pseudomonas aeruginosa.

Phage SN1 infects Sphaerotilus natans and Pseudomonas aeruginosa strains. Its genome consists of 61,858 bp (64.3% GC) and 89 genes, including 32 with predicted functions. SN1 genome is very similar to Pseudomonas phage M6, which contains hypermodified thymidines. Genome analyses revealed similar base-modifying genes as those found in M6.

A truncated anti-CRISPR protein prevents spacer acquisition but not interference.

CRISPR-Cas systems in prokaryotic cells provide an adaptive immunity against invading nucleic acids. For example, phage infection leads to addition of new immunity (spacer acquisition) and DNA cleavage (interference) in the bacterial model species Streptococcus thermophilus, which primarily relies on Cas9-containing CRISPR-Cas systems. Phages can counteract this defense system through mutations in the targeted protospacers or by encoding anti-CRISPR proteins (ACRs) that block Cas9 interference activity. Here, we show that S. thermophilus can block ACR-containing phages when the CRISPR immunity specifically targets the acr gene. This in turn selects for phage mutants carrying a deletion within the acr gene. Remarkably, a truncated acrIIA allele, found in a wild-type virulent streptococcal phage, does not block the interference activity of Cas9 but still prevents the acquisition of new immunities, thereby providing an example of an ACR specifically inhibiting spacer acquisition.

Complete Genome Sequences of 10 Lactococcal Skunavirus Phages Isolated from Cheddar Cheese Whey Samples in Canada.

We report the complete genome sequences of 10 virulent phages of the Skunavirus genus (Siphoviridae) that infect Lactococcus lactis strains used for cheddar cheese production in Canada. Their linear genomes range from 28,969 bp to 31,042 bp with GC contents of 34.1 to 35.1% and 55 to 60 predicted open reading frames (ORFs).

A Lactococcal Phage Protein Promotes Viral Propagation and Alters the Host Proteomic Response During Infection.

The lactococcal virulent phage p2 is a model for studying the Skunavirus genus, the most prevalent group of phages causing milk fermentation failures in cheese factories worldwide. This siphophage infects Lactococcus lactis MG1363, a model strain used to study Gram-positive lactic acid bacteria. The structural proteins of phage p2 have been thoroughly described, while most of its non-structural proteins remain uncharacterized. Here, we developed an integrative approach, making use of structural biology, genomics, physiology, and proteomics to provide insights into the function of ORF47, the most conserved non-structural protein of unknown function among the Skunavirus genus. This small phage protein, which is composed of three α-helices, was found to have a major impact on the bacterial proteome during phage infection and to significantly reduce the emergence of bacteriophage-insensitive mutants.

Source Tracking Based on Core Genome SNV and CRISPR Typing of Salmonella enterica Serovar Heidelberg Isolates Involved in Foodborne Outbreaks in Québec, 2012.

Whole-genome sequencing (WGS) is the method of choice for bacterial subtyping and it is rapidly replacing the more traditional methods such as pulsed-field gel electrophoresis (PFGE). Here we used the high-resolution core genome single nucleotide variant (cgSNV) typing method to characterize clinical and food from Salmonella enterica serovar Heidelberg isolates in the context of source attribution. Additionally, clustered regularly interspaced short palindromic repeats (CRISPR) analysis was included to further support this method. Our results revealed that cgSNV was highly discriminatory and separated the outbreak isolates into distinct clusters (0-4 SNVs). CRISPR analysis was also able to distinguish outbreak strains from epidemiologically unrelated isolates. Specifically, our data clearly demonstrated the strength of these two methods to determine the probable source(s) of a 2012 epidemiologically characterized outbreak of S. Heidelberg. Using molecular cut-off of 0-10 SNVs, the cgSNV analysis of 246 clinical and food isolates of S. Heidelberg collected in Québec, in the same year of the outbreak event, revealed that retail and abattoir chicken isolates likely represent an important source of human infection to S. Heidelberg. Interestingly, the isolates genetically related by cgSNV also harbored the same CRISPR as outbreak isolates and clusters. This indicates that CRISPR profiles can be useful as a complementary approach to determine source attribution in foodborne outbreaks. Use of the genomic analysis also allowed to identify a large number of cases that were missed by PFGE, indicating that most outbreaks are probably underestimated. Although epidemiological information must still support WGS-based results, cgSNV method is a highly discriminatory method for the resolution of outbreak events and the attribution of these events to their respective sources. CRISPR typing can serve as a complimentary tool to this analysis during source tracking.

DNA tandem repeats contribute to the genetic diversity of Brevibacterium aurantiacum phages.

This report presents the characterization of the first virulent phages infecting Brevibacterium aurantiacum, a bacterial species used during the manufacture of surface-ripened cheeses. These phages were also responsible for flavour and colour defects in surface-ripened cheeses. Sixteen phages (out of 62 isolates) were selected for genome sequencing and comparative analyses. These cos-type phages with a long non-contractile tail currently belong to the Siphoviridae family (Caudovirales order). Their genome sizes vary from 35,637 to 36,825 bp and, similar to their host, have a high GC content (~61%). Genes encoding for an immunity repressor, an excisionase and a truncated integrase were found, suggesting that these virulent phages may be derived from a prophage. Their genomic organization is highly conserved, with most of the diversity coming from the presence of long (198 bp) DNA tandem repeats (TRs) within an open reading frame coding for a protein of unknown function. We categorized these phages into seven genomic groups according to their number of TR, which ranged from two to eight. Moreover, we showed that TRs are widespread in phage genomes, found in more than 85% of the genomes available in public databases.

Novel Genus of Phages Infecting Streptococcus thermophilus: Genomic and Morphological Characterization.

Streptococcus thermophilus is a lactic acid bacterium commonly used for the manufacture of yogurt and specialty cheeses. Virulent phages represent a major risk for milk fermentation processes worldwide, as they can inactivate the added starter bacterial cells, leading to low-quality fermented dairy products. To date, four genetically distinct groups of phages infecting S. thermophilus have been described. Here, we describe a fifth group. Phages P738 and D4446 are virulent siphophages that infect a few industrial strains of S. thermophilus The genomes of phages P738 and D4446 were sequenced and found to contain 34,037 and 33,656 bp as well as 48 and 46 open reading frames, respectively. Comparative genomic analyses revealed that the two phages are closely related to each other but display very limited similarities to other S. thermophilus phages. In fact, these two novel S. thermophilus phages share similarities with streptococcal phages of nondairy origin, suggesting that they emerged recently in the dairy environment.IMPORTANCE Despite decades of research and adapted antiphage strategies such as CRISPR-Cas systems, virulent phages are still a persistent risk for the milk fermentation industry worldwide, as they can cause manufacturing failures and alter product quality. Phages P738 and D4446 are novel virulent phages that infect the food-grade Gram-positive bacterial species Streptococcus thermophilus These two related viruses represent a fifth group of S. thermophilus phages, as they are significantly distinct from other known S. thermophilus phages. Both phages share similarities with phages infecting nondairy streptococci, suggesting their recent emergence and probable coexistence in dairy environments. These findings highlight the necessity of phage surveillance programs as the phage population evolves in response to the application of antiphage strategies.

Author Correction: A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

A mutation in the methionine aminopeptidase gene provides phage resistance in Streptococcus thermophilus.

Streptococcus thermophilus is a lactic acid bacterium widely used by the dairy industry for the manufacture of yogurt and specialty cheeses. It is also a Gram-positive bacterial model to study phage-host interactions. CRISPR-Cas systems are one of the most prevalent phage resistance mechanisms in S. thermophilus. Little information is available about other host factors involved in phage replication in this food-grade streptococcal species. We used the model strain S. thermophilus SMQ-301 and its virulent phage DT1, harboring the anti-CRISPR protein AcrIIA6, to show that a host gene coding for a methionine aminopeptidase (metAP) is necessary for phage DT1 to complete its lytic cycle. A single mutation in metAP provides S. thermophilus SMQ-301 with strong resistance against phage DT1. The mutation impedes a late step of the lytic cycle since phage adsorption, DNA replication, and protein expression were not affected. When the mutated strain was complemented with the wild-type version of the gene, the phage sensitivity phenotype was restored. When this mutation was introduced into other S. thermophilus strains it provided resistance against cos-type (Sfi21dt1virus genus) phages but replication of pac-type (Sfi11virus genus) phages was not affected. The mutation in the gene coding for the MetAP induces amino acid change in a catalytic domain conserved across many bacterial species. Introducing the same mutation in Streptococcus mutans also provided a phage resistance phenotype, suggesting the wide-ranging importance of the host methionine aminopeptidase in phage replication.

Diversity and Host Specificity Revealed by Biological Characterization and Whole Genome Sequencing of Bacteriophages Infecting Salmonella enterica.

Phages infecting members of the opportunistic human pathogen, Salmonella enterica, are widespread in natural environments and offer a potential source of agents that could be used for controlling populations of this bacterium; yet, relatively little is known about these phages. Here we describe the isolation and characterization of 45 phages of Salmonella enterica from disparate geographic locations within British Columbia, Canada. Host-range profiling revealed host-specific patterns of susceptibility and resistance, with several phages identified that have a broad-host range (i.e., able to lyse >40% of bacterial hosts tested). One phage in particular, SE13, is able to lyse 51 out of the 61 Salmonella strains tested. Comparative genomic analyses also revealed an abundance of sequence diversity in the sequenced phages. Alignment of the genomes grouped the phages into 12 clusters with three singletons. Phages within certain clusters exhibited extraordinarily high genome homology (>98% nucleotide identity), yet between clusters, genomes exhibited a span of diversity (<50% nucleotide identity). Alignment of the major capsid protein also supported the clustering pattern observed with alignment of the whole genomes. We further observed associations between genomic relatedness and the site of isolation, as well as genetic elements related to DNA metabolism and host virulence. Our data support the knowledge framework for phage diversity and phage-host interactions that are required for developing phage-based applications for various sectors, including biocontrol, detection and typing.

Variability in the durability of CRISPR-Cas immunity.

The durability of host resistance is challenged by the ability of pathogens to escape the defence of their hosts. Understanding the variability in the durability of host resistance is of paramount importance for designing more effective control strategies against infectious diseases. Here, we study the durability of various clustered regularly interspaced short palindromic repeats-Cas (CRISPR-Cas) alleles of the bacteria Streptococcus thermophilus against lytic phages. We found substantial variability in durability among different resistant bacteria. Since the escape of the phage is driven by a mutation in the phage sequence targeted by CRISPR-Cas, we explored the fitness costs associated with these escape mutations. We found that, on average, escape mutations decrease the fitness of the phage. Yet, the magnitude of this fitness cost does not predict the durability of CRISPR-Cas immunity. We contend that this variability in the durability of resistance may be because of variations in phage mutation rate or in the proportion of lethal mutations across the phage genome. These results have important implications on the coevolutionary dynamics between bacteria and phages and for the optimal deployment of resistance strategies against pathogens and pests. Understanding the durability of CRISPR-Cas immunity may also help develop more effective gene-drive strategies based on CRISPR-Cas9 technology. This article is part of a discussion meeting issue 'The ecology and evolution of prokaryotic CRISPR-Cas adaptive immune systems'.

Complete Genome Sequence of Escherichia coli Siphophage BRET.

The lytic Escherichia coli siphophage BRET was isolated from a chicken obtained at a local market in Abidjan, Côte d'Ivoire. Its linear genome sequence consists of 59,550 bp (43.4% GC content) and contains 88 predicted genes, including 4 involved in archaeosine biosynthesis. Phage BRET is related (95% nucleotide identity) to Enterobacteria phage JenK1.

Lactococcus lactis type III-A CRISPR-Cas system cleaves bacteriophage RNA.

CRISPR-Cas defends microbial cells against invading nucleic acids including viral genomes. Recent studies have shown that type III-A CRISPR-Cas systems target both RNA and DNA in a transcription-dependent manner. We previously found a type III-A system on a conjugative plasmid in Lactococcus lactis which provided resistance against virulent phages of the Siphoviridae family. Its naturally occurring spacers are oriented to generate crRNAs complementary to target phage mRNA, suggesting transcription-dependent targeting. Here, we show that only constructs whose spacers produce crRNAs complementary to the phage mRNA confer phage resistance in L. lactis. In vivo nucleic acid cleavage assays showed that cleavage of phage dsDNA genome was not detected within phage-infected L. lactis cells. On the other hand, Northern blots indicated that the lactococcal CRISPR-Cas cleaves phage mRNA in vivo. These results cannot exclude that single-stranded phage DNA is not being targeted, but phage DNA replication has been shown to be impaired.

Characterization of CRISPR-Cas systems in the Ralstonia solanacearum species complex.

Clustered regularly interspaced short palindromic repeats (CRISPRs) are composed of an array of short DNA repeat sequences separated by unique spacer sequences that are flanked by associated (Cas) genes. CRISPR-Cas systems are found in the genomes of several microbes and can act as an adaptive immune mechanism against invading foreign nucleic acids, such as phage genomes. Here, we studied the CRISPR-Cas systems in plant-pathogenic bacteria of the Ralstonia solanacearum species complex (RSSC). A CRISPR-Cas system was found in 31% of RSSC genomes present in public databases. Specifically, CRISPR-Cas types I-E and II-C were found, with I-E being the most common. The presence of the same CRISPR-Cas types in distinct Ralstonia phylotypes and species suggests the acquisition of the system by a common ancestor before Ralstonia species segregation. In addition, a Cas1 phylogeny (I-E type) showed a perfect geographical segregation of phylotypes, supporting an ancient acquisition. Ralstoniasolanacearum strains CFBP2957 and K60T were challenged with a virulent phage, and the CRISPR arrays of bacteriophage-insensitive mutants (BIMs) were analysed. No new spacer acquisition was detected in the analysed BIMs. The functionality of the CRISPR-Cas interference step was also tested in R. solanacearum CFBP2957 using a spacer-protospacer adjacent motif (PAM) delivery system, and no resistance was observed against phage phiAP1. Our results show that the CRISPR-Cas system in R. solanacearum CFBP2957 is not its primary antiviral strategy.

Characterization of the Escherichia coli Virulent Myophage ST32.

The virulent phage ST32 that infects the Escherichiacoli strain ST130 was isolated from a wastewater sample in China and analyzed. Morphological observations showed that phage ST32 belongs to the Myoviridae family, as it has an icosahedral capsid and long contractile tail. Host range analysis showed that it exhibits a broad range of hosts including non-pathogenic and pathogenic E. coli strains. Interestingly, phage ST32 had a much larger burst size when amplified at 20 °C as compared to 30 °C or 37 °C. Its double-stranded DNA genome was sequenced and found to contain 53,092 bp with a GC content of 44.14%. Seventy-nine open reading frames (ORFs) were identified and annotated as well as a tRNA-Arg. Only nineteen ORFs were assigned putative functions. A phylogenetic tree using the large terminase subunit revealed a close relatedness with four unclassified Myoviridae phages. A comparative genomic analysis of these phages showed that the Enterobacteria phage phiEcoM-GJ1 is the closest relative to ST32 and shares the same new branch in the phylogenetic tree. Still, these two phages share only 47 of 79 ORFs with more than 90% identity. Phage ST32 has unique characteristics that make it a potential biological control agent under specific conditions.